Anomalous double peak structure in Nb/Ni superconductor/ferromagnet tunneling 
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We have experimentally investigated the density of states (DOS) in Nb/Ni (S/F) bilayers as a 
function of Ni thickness, dp. Our thinnest samples show the usual DOS peak at ±Ao, whereas 
intermediate-thickness samples have an anomalous "double-peak" structure. For thicker sam- 
ples (d_F > 3.5 nm), we see an "inverted" DOS which has previously only been reported in 
superconductor/ weak- ferromagnet structures. We analyze the data using the self-consistent non- 
linear Usadel equation and find that we are able to quantitatively fit the features at ±Ao if we 
include a large amount of spin-orbit scattering in the model. Interestingly, we are unable to re- 
produce the sub-gap structure through the addition of any parameter(s). Therefore, the observed 
anomalous sub-gap structure represents new physics beyond that contained in the present Usadel 
theory. 



The co-existence of superconductivity and ferromag- 
netism was first proposed by Fulde and Ferrell [l[ and 
Larkin and Ovchinnikov [2] more than forty years ago. 
While some unusual materials have since been found with 
both superconducting and ferromagnetic transitions {e.g. 
ErRli4B4 [3|), much recent interest has focused on con- 
ventional superconductor/ferromagnet (S/F) proximity 
effect multi-layer systems. A wide variety of unusual phe- 
nomena has been proposed for these systems including 
oscillating critical temperatures, T c 0], 7r-state Josesph- 
son junctions and long-ranged odd- frequency triplet 
superconductivity 

Qualitative evidence for the first two of these effects 
is convincing, but definitive quantitative agreement with 
theory has been problematic. The evidence for triplet 
superconductivity is less certain, although a recent re- 
port by Keizer et al provides tantalizing evidence for 
such an effect. One reason for the difficulty in achieving 
quantitative agreement with theory is the proliferation of 
physical effects that now have been incorporated into the 
theory, leading to a concomitant proliferation of fitting 
parameters, which makes discriminating fits to limited 
data sets nearly impossible. 

In order to obtain more discriminating data sets and 
to further explore the SF proximity effect in the case of 
strong ferromagnets, we have undertaken superconduct- 
ing tunneling densities of state (DOS) measurements on 
Nb/Ni thin-film bilayers. By varying the Ni thickness, 
dp, we can track the spatial evolution of the behavior of 
the Cooper pairs diffusing into the ferromagnet. This ap- 
proach gives us much more information per sample (the 
entire DOS spectrum) than T c or J c measurements, and is 
less sensitive to variations in boundary parameters. Ana- 
lyzing these results with the most complete forms of the 
Usadel theory available has allowed us to discriminate 



critically for the first time the relative importance of the 
various physical effects now incorporated into the theory. 
We find that by far the most important parameter be- 
yond the exchange field, E ex , is the degree of spin-orbit 
scattering (first suggested by Demler et al [1]). In addi- 
tion, we find an anomalous double-peak structure in the 
DOS that has not been reported previously and that we 
have been unable to account for theoretically. 

We use planar tunnel junctions of the form normal- 
insulator-ferromagnet-superconductor. A schematic of 
our sample geometry is shown in the inset of Fig. [2] 
The deposition of our samples and characterization of the 
tunnel junctions has been documented elsewhere (jl. [Toj] . 
In brief, the various layers are sputtered and patterned 
with stencil masks in situ in a DC magnetron sputtering 
chamber without breaking vacuum. The AI2O3 tunnel 
barriers are formed by oxidizing the Al underlayer, and 
are canonical in their behavior, except for the zero-bias 
anomalies commonly observed in tunnel junctions incor- 
porating magnetic materials; we discuss these below. We 
use a Co6oFe4o backing of the Al electrode so as to reduce 
its critical temperature below our lowest measurement 
temperature, which ensures normal/superconductor tun- 
neling. We vary oIf from nm to 5 nm in 0.5 nm incre- 
ments. Ni has a Curie Temperature of roughly 600 K, 
and should be ferromagnetic for film thicknesses greater 
than two atomic layers (5 A) [III, |T2|. We have taken 
care to ensure that our Ni films are ferromagnetic by 
measuring the superconducting resistive transition of the 
bilayer in a parallel magnetic field: we detect hysteretic 
signals confirming the presence of magnetism in all our 
films with dp > 1.0 nm. 

Tunneling measurements were performed at 0.28 K us- 
ing standard lock-in techniques. Samples were measured 
in zero magnetic field and also in a perpendicular field 
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greater than H C 2. In all of our samples, the normal-state 
conductance is a V-shaped curve in which the zero-bias 
conductance (ZBC) is roughly 1% — 2% lower than the 
conductance at 5 mV, which is typical of magnetic tun- 
nel junctions (l3j . Since there are no systematic trends in 
the size of the background, we do not consider it relevant 
to the superconducting DOS. As discussed in more detail 
m Ref. 1 101, we remove this background conductance in 
our data analysis by dividing the zero-field conductance 
by the high-field conductance, thus isolating the super- 
conducting DOS. Using this procedure, we find that the 
resultant DOS satisfies the sum rule on the total density 
of states, except in our thickest samples, i.e. those with 
the smallest conductance variations. 



The results of this normalization procedure for all mea- 
sured dp are shown in Fig. [T] with dp — (tunneling 
into pure Nb) at the bottom and dp = 5.0 nm at the 
top, where the conductance scale has been magnified by 
roughly a factor of 1000 in comparison. The dp = curve 
has two clear BCS coherence peaks at ±Ao = 1.3 meV; 
above that, we see that the addition of just 1.0 nm of Ni 
increases the ZBC significantly. In addition to the gen- 
eral trend of decreasing feature size as a function of dp, 
for cLf = 1-5 nm we see a striking new feature: two co- 
herence peaks on either side of zero-bias. Continuing up 
Fig. [TJ we see this "double-peak" structure of the DOS 
for dp = 1.5 — 3.0 nm, with the interior peaks moving to 
lower voltages as dp increases. At dp = 3.5 nm, though, 
we see a qualitatively different DOS: the coherence peaks 
at V = ±Ao have "inverted" and are now conductance 
minima, while a narrow interior gap remains. This in- 
verted DOS was seen in the first reported measurement 
of the DOS in S/F bilayers 14[, yet has not been seen in 
other reports [lO, EH ■ The outer peak remains inverted 
at dp = 4.0 nm, but by dp = 4.5, 5.0 nm we have re- 
turned to the non-inverted DOS with the addition of a 
very narrow peak at zero- bias. These measurements of 
thicker samples (dp > 4.5 nm) must be taken with a 
grain of salt as the feature size is quite small. It is our 
opinion that this zero-bias peak, which - like all other 
features of the DOS - does not change in a small par- 
allel or perpendicular magnetic field is due to the steep 
voltage dependence of the background conductance and 
is therefore a by-product of our normalization procedure. 



We model our system in the dirty limit, which is ap- 
plicable only when the elastic scattering time, r e , is 
shorter than all other relevant time-scales. In the F- 
layer, H/r e « 400 meV (from resistivity data), which is 
much greater than either the estimated exchange field, 
E ex = 78 meV (estimated from the Curie Tempera- 
ture [HI]), or the superconducting gap, A = 1.3 meV. In 
this limit, the superconducting order parameter should 
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FIG. 1: Normalized conductances taken at 0.28 K for various 
Ni thickness indicated inside. From the bottom to the top 
plot the vertical scale is successively amplified. The curves 
are shifted for clarity. The vertical lines show that the outside 
peaks energy remains unchanged from 0.5 to 4.5 nm. 



obey the Usadel equation [17|, [1 



~ ~2 — !h? + (~ iu ± lEex + 2Tz cos T(i) ) sin T(i)+ 

T x sin(# T + t ) ± r so sin(# T - 0;) = Acos0 TU) , (1) 

where 0-m) corresponds to the up (down) electron den- 
sity, D is the diffusion constant, ui is the energy, E ex is the 
exchange field, T z and Tx are the strength of the mag- 
netic scattering parallel and perpendicular to the quan- 
tization axis, and T$o is the strength of the spin-orbit 
scattering. This equation is valid in the S-layer when 
all magnetic terms (E ex , Tz, Tx) are zero and in the F 
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layer when the gap, A, is zero. A must also obey the 
self-consistent equation, which at zero temperature is 



A(a;) = A 



dw(l/2)Im [sin ( 



HI 



(2) 



where uijj is the Debye frequency, and A is the coupling 
constant; this can be solved iteratively. 

The S-layer and F-layer occupy the space —ds < x < 
and < x < d F , respectively. We supplement our equa- 
tion with the usual (non-magnetic) boundary conditions: 
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where £, N = y/(hD F /2A), ? s = y/(hD 3 /2A), lB = 
RbA/ pf£n, and 7 = _£s£,s / Pf£,n ■ This notation is mod- 
eled after Refs. [la , Il9| . but with minor differences for 
ease of calculation. Finally, the total DOS measured by 
our junctions, N(u>), is 



N(u>) = (l/2)Re [cos Of (w) + cos9f (w)] x=d 



(5) 



If we numerically solve the Usadel equation with 
E ex 3> A and no scattering, we find the characteris- 
tic decay and oscillation of the DOS as a function of 
d F [Uli but as dp increases, it goes from inverted to 
normal, to inverted, etc. In other words, we expect the 
inverted — not the normal — DOS to be seen for thin 
(0.2 nm < dp < 5 nm) F layers. Cottet et al 2l| were 
the first to comment on this unusual prediction — which 
has not been seen in any tunneling study of S/F systems 
- and suggested that spin-dependent interfacial phase 
shifts (SDIPS) could resolve this disagreement. 

We find that we can qualitatively account for the ob- 
served behavior through the addition of either Tz, T$0: 
or SDIPS — the quantitative fits discussed below favor 
Tso a s the dominant scattering term. Equally impor- 
tant, as noted before, we find that once the Usadel equa- 
tion in the F-layer is essentially linear (dp f» £p) the 
only effect of increasing dp is to scale the DOS at the 
interface [13, H3] • This means that none of these param- 
eters will produce any sub-gap structure - nor, for that 
matter, will any other parameter in Eq. [T] 25] 

Thus, we are motivated to split empirically our DOS 
into two parts: a DOS with a single peak (or inverted 
peak) at ±Ao, which we call the outer gap, and what- 
ever remains, which we call the sub gap. Since the 
dp = 1.0 nm curve only has peaks at ±Ao we can 
use it as a template to isolate the outer-gap features 
of the other curves. We want to break down the other 
curves into a sum of two curves: one that is a scaled 
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FIG. 2: (a) Normalized conductance of the dp = 2.0 nm 
sample, split into outer-gap and sub-gap contributions. The 
outer-gap contribution (dashed line), a scaled version of the 
dp = 1.0 nm curve and the sub-gap contribution (dash-dot 
line) sum to the total conductance (solid line) . (b) Normalized 
conductance of the d F — 1.5 nm sample taken in a perpen- 
dicular magnetic field, split into outer and sub-gap contribu- 
tions. The sub-gap peak shifts in energy and changes ampli- 
tude slightly, but it does not broaden in an applied field, like 
the outer-gap. 



(in the conductance-axis) copy of the template and an- 
other which will contain all of the sub-gap features. In 
order to determine the size of the different contributions, 
we adjust the scaling of the template such that the re- 
maining sub-gap contribution is as smooth as possible at 
±Aq. Figure[2] (a) shows the results of this process on 
the dp = 2.0 nm sample. Qualitatively, the anomalous 
contribution for 1.5 < dp < 3.5 nm looks like a super- 
conducting gap whose width decreases with increasing dp 
and which stays un-inverted even when outer gap inverts. 

Once we have isolated the outer-gap contribution for 
all dp, we can analyze it quantitatively with the Usadel 
equation. As noted above, a combination of E ex , SDIPS, 
Tz, and Tso could be used, but we find that the best 
and simplest fits require only E ex and Tso ■ Figure [3] 
shows the absolute value of 1— N(0), the ZBC, versus dp. 
The circles represent the ZBC of the outer-gap contribu- 
tion and the triangles represent the sub-gap contribution. 
The line is the calculated ZBC from the Usadel equa- 
tion with parameters E ex = 63 meV, Tso — 46 meV, 
7s = 0.54, and 7 = 0.52, which fits the normal ZBC rea- 
sonably well [2^1 . Thus, we conclude that the outer gap is 
well-described by the Usadel equation with an exchange 
field and spin-orbit scattering. 

We measured the d F = 1.5 nm sample in a perpen- 
dicular magnetic field from B = to 3000 G in 500 G 
increments. To isolate the outer gap, we measured the 
template curve {dp =1.0 nm) at the same fields and used 
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FIG. 3: Abs(l - AT(0)), the zero-bias conductance (ZBC), 
split into outer-gap (circles) and sub-gap (triangles) contribu- 
tions. Open (filled) markers represent a ZBC less (greater) 
than unity. The outer-gap contribution is fitted to Eq.[T]with 
Eex = 63 meV and Vso = 46 meV (thin line). 

those as templates for each applied field. Figure [2] (b) 
shows the results of this process: the outer curves are 
scaled versions of the template while the inner curves are 
the remaining sub-gap feature. The relative weight of the 
template needed to isolate each sub-gap curve is roughly 
constant, but the resulting sub-gap curves decrease in 
size with increasing field. Further, the shape of the sub- 
gap curve remains virtually unchanged in field, while the 
outer-gap peaks display significant broadening. 

The existence of sub-gap structure alone indicates 
physics beyond the standard Usadel treatment. It is in- 
teresting to speculate what that might be. One intriguing 
possibility is that it is related to triplet superconductiv- 
ity, perhaps in combination with some spin-flip scattering 
at the SF interface (see Bergeret et al [23| and references 
therein) . While it is clearly premature to make this case 
confidently, we note that the relative insensitivity of the 
shape of the sub-gap peak to magnetic field and the ap- 
parent slower decay of its associated ZBC as a function 
of d,F are qualitatively consistent with such a possibil- 
ity. Additional theory will be required to fully assess this 
possibility. 

In summary, we have measured the DOS of Nb/Ni bi- 
layers as a function of Ni thickness. In addition to tun- 
neling features which are well-described by the Usadel 
equation with a strong exchange field and spin-orbit cou- 
pling, we have also discovered a robust sub-gap structure 
which cannot be explained by the conventional theory. 
By isolating this sub-gap contribution, we have shown 
that it behaves in a qualitatively different way from the 
Usadel contribution in a perpendicular magnetic field, 
which leads us to inquire whether it is the result of m = 1 



triplet correlations in the bilayers. 
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